Spatially heterogeneous dynamics in granular compaction 
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We prove the emergence of spatially correlated dynamics in slowly compacting dense granular 
media by analyzing analytically and numerically multi-point correlation functions in a simple particle 
model characterized by slow non-equilibrium dynamics. We show that the logarithmically slow 
dynamics at large times is accompanied by spatially extended dynamic structures that resemble the 
ones observed in glass-forming liquids and dense colloidal suspensions. This suggests that dynamic 
heterogeneity is another key common feature present in very different jamming materials. 
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When a pile of grains is gently shaken, its volume frac- 
tion increases so slowly that the process hardly becomes 
stationary on experimental timescales E10- This is rem- 
iniscent of the slow relaxation observed in glass-formers, 
as noted long ago Q . From a more fundamental point of 
view, it is tempting to build upon analogies and suggest 
that granular media, glasses, and other jamming systems 
can be described by common theoretical approaches Q. 
In recent years, several aspects of glasses and granular 
media have been studied with similar approaches. Static 
structures have been studied to understand the relevance 
of the network of force chains between grains or atoms 
to dynamics of the jammed state Also, since both 
grains and glasses undergo non-equilibrium 'glassy' dy- 
namics, the idea that an effective thermodynamics can 
be used has received considerable attention [1 . 

In this work we also transfer knowledge from one field 
to the other and show that the glassy dynamics of gran- 
ular media is characterized by the appearance of spatio- 
temporal structures similar to the ones described as dy- 
namic heterogeneity in glass-formers and dense col- 
loidal suspensions [!j . Dynamic heterogeneity is believed 
to play a crucial role in the glass formation, and forms 
the core of recent theoretical descriptions [lfj. Physi- 
cally, it stems from the existence of spatial correlations in 
the local dynamics that extend beyond the ones revealed 
by static pair correlations. To study dynamic hetero- 
geneity, correlators that probe more than two points in 
space and time have to be considered d 0, ^| . These 
spatial fluctuations have never been studied in models 
or experiments on granular compaction, although caged 
particle dynamics was recently studied in a sheared sys- 
tem |l2j . Here we prove the emergence of large dynamic 
lengthscales in a particle model which is a variant of the 
parking lot model, introduced and studied in detail in the 
context of granular compaction 
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take advantage of its relative simplicity to compute ana- 
lytically multi-point correlations studied in glass-formers 
and confirm our results by numerical simulations. 

We consider a variant of the parking lot model E3| in- 
troduced in Ref. [1J|. The model is a one-dimensional 



process in which hard blocks of unit size are first irre- 
versibly deposited at random positions on a line of linear 
size L until no place is available for more depositions. 
Timescales are counted from the end of the deposition 
process which corresponds to t w = 0. In a second step, 
particles are allowed to diffuse with a coefficient of dif- 
fusion D. The last dynamical rule defining the model 
consists of deposition events. When particles have dif- 
fused in such a way that a void of unit size opens, the 
void is instantaneously filled by a new particle and the 
density, p(t w ) — N(t w )/L, increases by 1/L, iV(£ w ) being 
the total number of particles present at time t w . These 
rules lead to slow dynamics, because the larger the den- 
sity the longer it takes to open a void. A more general 
version of this model includes evaporation 0, 0. It 
displays the generic features observed in granular com- 
paction: logarithmic increase of the density [l^l lldl llfij . 
aging Ealrjjj non-Gaussian density fluctuations |lj| .ef- 
fective temperatures and link with Edwards entropy [15( , 
hysteresis effects Q. The time evolution of the density 
also describes well the experiments, 
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where poo, ao, a\, r and to are fitting parameters. For the 
present model without evaporation Eq. Q holds at large 
i w with poo = a\ = t = 1, ao = to = 0. Let us remark 
here that such logarithmic compaction has also been ob- 
tain from different models with depositiomevaporation 
and diffusion, like those studied in Ref. 18] . Here we 
focus on local dynamic quantities and their spatial cor- 
relations through both analytic calculations and Monte 
Carlo simulations. The results of the simulations have 
been obtained by averaging over 2.10 4 independent his- 
tories with L = 250. Time is measured in units of the 
diffusion constant, D, which is set to unity. 

In Fig. we illustrate the results reported in this paper, 
namely the emergence of spatio-temporal correlations in 
dynamic trajectories of compacting granular systems. In 
this figure, black (white) denote particles that are more 
(less) mobile than the average on a timescale of the order 
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FIG. 1: Space-time pattern of dynamic fluctuations: par- 
ticles are represented in black (white) when moving slower 
(faster) than average. Local dynamics is measured via the 
self-intermediate scattering function, Eq. 0, with k = it, at 
fixed i-iw = 2.10 4 . 



of the structural relaxation time. The time extension of 
the trajectory is 10 5 , the spatial extension L = 500, and 
the final density p(l0 5 ) ~ 0.91. The emergence of large- 
scale correlations in the local dynamics is evident, since 
domains that extend spatially over a hundred particles 
can be observed, despite the absence of any such large 
static correlations between particles. Spatial clustering 
of fast and slow reg ions is the hallmark of dynamic het- 
erogeneity |R|jj jlCllTl| . A closer look at the left part of 
the figure shows that the dynamic lengthscale is visibly 
smaller at shorter times, suggesting that spatial correla- 
tions are larger when the dynamics becomes slower, and 
possibly diverge in the jamming limit of a full system. 

How did we build Fig. \I$ As in supercooled liquids, 
information on local dynamics is accessed by following 
the dynamics of individual particles. Consider two times 
t w and t > t w 3> 1. The distribution of particle displace- 
ments is the self-part of the van-Hove function, 
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where Sri(t,t w ) = rj(i) — ri(t w ) is the displacement of 
particle i between times t w and t and only particles 
present at time t w are summed over. Its Fourier trans- 
form is the self-intermediate scattering function 

F s (k,t,t w ) = — — (cos(kSn(t,t w ))). (3) 

Aging will manifest itself through an explicit dependence 
of F s (k,t,t w ) on its two time arguments In the fol- 



lowing we focus on the value k = 2ir which corresponds 
to displacements of the order of the particle size. 

Since dynamic heterogeneity relates to spatial fluctu- 
ations about the averaged two-time dynamics, we have 
shown in Fig. ^ spatial fluctuations about F s (k,t,t w ). 
We colored black (white) those particles for which 
5Fi(k,t,t w ) = cos (kSri(t, t w )) — F s (k,t,t w ) is negative 
(positive), i.e. those particles that move more (less) than 
average in a particular realization. 

To further quantify how particle displacements are spa- 
tially correlated we consider the structure factor of the 
dynamic heterogeneity shown in Fig. ^ 
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where fk(t, t w ) — (5Fi(k, t, t w ) 2 } normalizes the structure 
factor and some obvious time and wavevector dependen- 
cies have been removed, for clarity. By definition, the 
structure factor is built from two-time two-point quan- 
tities which is the minimum requirement to probe the 
spatio-temporal patterns of Fig. ^ The small q limit in 
Eq. J3J defines a dynamic "four-point" susceptibility, 



Xk(t, t w ) = S k (q = 0, t, t w ), 



(5) 



which can be rewritten as the variance of the fluctua- 
tions of the spatially averaged two-time dynamics. Phys- 
ically, dynamic fluctuations increase when the number of 
independently relaxing objects decreases Normal- 
izations ensure that \k is finite in the thermodynamic 
limit, except at a dynamic critical point [Tsj. 

The quantities introduced above can be estimated an- 
alytically as follows. Between two deposition events, the 
total free volume is L(l — p) and thus the slow dynamics 
can be shown |20j | to be equivalent to a continuum version 
of asymmetric exclusion process at density p — p/(l — p), 
i.e. involving caged dynamics. This process is then 
mapped to a fluctuating interface model 0, H2- 
Labeling particles in the order of their position from 
the origin we introduce an interface position y(x,t) for 
the ith particle encountered on the lattice, such that 
y(x,t) = ri(t), x replacing i as a label for the parti- 
cle. We then write y(x,t) = h(x,t) + x/p, were h(x,t) 
is the deviation from the position of the particle from an 
uniform configuration. The connection to symmetric ex- 
clusion process allows one to show [13, El El that h(x, t) 
follows an Edwards- Wilkinson dynamics |23j, 
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where a — p 2 ~ ln(£ w ) 2 and 77 is a Gaussian white noise 
with zero mean and variance (r)(x, t)r)(x' , if)) = 28(x — 
x')S(t — t'). From © the equilibrium measure reads 



P cq [h(x)} = exp (-^- J dx [Vh{ 
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which shows that deviations from the averaged position 
of a particle arise with an elastic penalty in the interface 
representation. This elastic behavior makes the following 
calculation similar to the evaluation of the elastic contri- 
bution to the dynamical susceptibility JSJ in supercooled 
liquids performed in Ref. |24j . 

Equation © is first solved in the Fourier space, 
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where we have defined the time difference r — t — i w . A 
crucial approximation is made here since a in Eq. JSJ has 
in fact a logarithmic time dependence. This amounts to 
neglecting the effect of further deposition events on the 
particles already present at t w . Since deposition is such a 
rare event at large times this approximation should cap- 
ture the evolution of the local dynamics, as our numerical 
simulations shall confirm. Using JSJ, one easily gets 



F s (k,t,t w ) = exp -fc 2 
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FIG. 2: Comparison of the numerical and analytical evalu- 
ations of Xk(~t,tw) as a function of r = t — t„ for k = 2n 
and t w = 6250, 12500, 25000, 50000, 100000 (from bottom to 
top). 



where g(x) = f™ x dq e^ x {\ - e"? 2 / 2 )/^, so that g(0) = 
(27T) 1 / 2 . Equation © is a classic result for the ordering 
process of one-dimensional random walkers (2^, [2(| . It 
shows that F s (k,t,t w ) displays aging behavior and scales 
with T/Tfc(i w ), with a logarithmically increasing relax- 
ation timescale Tk(t w ) ~ a(t w )/k 4 . Particles are there- 
fore sub-diffusing, (Srf(t,t w )) ~ y/r/a(t w ), and the van 
Hove function (J2J) is a Gaussian. Gaussiannity is a con- 
sequence of the diffusive nature of the particle motion 
at short times, which explains the discrepancy with the 
non-Gaussian distributions of displacements that have re- 
cently been measured experimentally in a bidimensional 
geometry [12]. 

The dynamic susceptibility (J5J can also be computed, 
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where T(x) is a scaling function defined from the first 
line of Eq. (|10(l . Careful analysis of Eq. (|10fl shows that 
Xk(t, t w ) goes from zero at r = to the asymptotic value 



Xk(t 



,*w) = a(t w )/(2fc 2 | 5 '(0)|), 



via a maximum 



Xl ~ a(t w )/k 2 when r - r fe (i w ). 

In Fig. [21 we show the dynamic susceptibility as a func- 
tion of time separation r for various ages t w obtained 
from direct numerical simulations. The lines through the 
data are from our analytical estimation, Ea. (|10J) . which 
constitutes clearly an excellent quantitative representa- 
tion of the data when t w becomes large. This supports 
the approximation made above that the dynamics at 



large times is slow enough that deposition events have 
little influence on local dynamics. 

The curves in Fig. [21 strikingly resemble the four- 
point susceptibilities discussed for supercooled liquids ap- 
proachin g th eir glass transition || and simple coarsening 
systems [ll|. As for those systems, we conclude that 
dynamics is maximally heterogeneous when observed at 
time separations close to the relaxation timescale, itself 
dependent of the age of the sample. To our knowledge 
no experimental determination of the dynamic suscep- 
tibility Xk(t,t w ) has been reported for granular media, 
although the experimental set-up described in Ref. 0] 
would probably allow its determination. 

The dynamics susceptibility Xk(t, t w ) measures the vol- 
ume integral of a spatial correlator. Therefore, an in- 
creasing susceptibility directly suggests the existence of 
a growing dynamic correlation length. This is most di- 
rectly seen in the Fourier space when the wavevector de- 
pendence of the dynamic structure factor Q is consid- 
ered. We have obtained an analytical form for Sk(q, t, t w ) 
but it is too lengthy to be reported here [2(J ■ In Fig. we 
show the wavector dependence of Sk(q, t, t w ) for various 
t w at time separations corresponding to the maximum of 
the dynamic susceptibility, obtained in numerical simu- 
lations. We also show the analytical results. As for the 
dynamic susceptibility the agreement between analytical 
and numerical results is excellent. 

At fixed t w the structure factor is characterized by a 
plateau at small q whose height is given by Xk(t w ). When 
q increases Sk(q, t, t w ) leaves the plateau and decreases to 
at large q. When t w is increased the plateau becomes 
higher and it ends at a smaller wavevector but its overall 
shape is unchanged. This implies a dynamic scale invari- 
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FIG. 3: Inset: Simulated dynamic structure factor Q as a 
function of q for parameters as in Fig. [5] and for r = Th(t v ). 
Main: Dynamic scale invariance, Eq. is revealed by 

rescaling space by £fc and Sk by x£- The full line is the an- 
alytical result using the mapping to the Edwards- Wilkinson 
interface model. 



ance such as found in glass- formers [ToL Il9j |: rescaling 
times by Tk(t w ) and space by £fc(t w ) ~ Oi{t w )/k 2 makes 
trajectories statistically equivalent. Formally this means 
that the following scaling law is obeyed, 



S k {q,T k ,t w ) = x* k (t w )Q k (g^fe(tw)) • 



(11) 



Again we note that data in Fig. [3] strikingly resemble 
dynamic structure factors measured both in realistic su- 
percooled liquids lal and coarse-grained models for the 
glass transition jic| . 

A major result of the above analysis is the existence of 
a dynamic length scale, £fc(i w ), which grows logarithmi- 
cally with time when compaction proceeds, and therefore 
diverges when the systems jams. A diverging lengthscale 
provides support to the temporal renormalization group 
argument developed in Ref. 0], but we see no obvi- 
ous connection between the dynamic criticality described 
here and the various power law scalings observed in static 
systems approaching jamming from above |27| . Physi- 
cally, collective rearrangements of particles are needed to 
create a void of unit size, the denser the system the more 
cooperative the dynamics. A naive determination of £ 
would rely on a free volume argument |28j . The mean 
free volume available to particle is so that the num- 
ber of particles required to have a fluctuation of the free 
volume equal to unity is Af ~ p 2 ~ £,{t w )- This simple 
physical argument underlies the cooperative nature of the 
dynamics which is more formally captured by four-point 
correlation functions, Eqs. Q and 

The idea that the size of collective motions increases 
when dynamics becomes slow is certainly not new. Multi- 
point correlation functions have now been measured in 



very different materials with similar qualitative results. 
That large dynamic lengthscales control glassy dynam- 
ics suggests the possibility that few universality classes 
underly and possibly unify the dynamical behavior of a 
much wider diversity of jamming materials. 

We thank S.N. Majumdar for enlightening discussions. 
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